Fast, high fidelity information transmission through spin chain quantum wires 
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Spin chains have been proposed as quantum wires for information transfer in solid state quantum 
architectures. We show that huge gains in both transfer speed and fidelity are possible using a 
minimalist control approach that relies only a single, local, on-off switch actuator. Effective switching 
time sequences can be determined using optimization techniques for both ideal and disordered chains. 
Simulations suggest that effective optimization is possible even in the absence of accurate models. 



The dynamical evolution of a chain of coupled spins is 
well-suited to the job of a data bus, transferring quan- 
tum information over the short distances between quan- 
tum registers. This is particularly desirable in architec- 
tures where internal quantum communication is required 
but the use of photons is impractical [TJ [SJ. A practi- 
cal scalable data-bus should achieve fast, high fidelity of 
information transfer without the need for intricate ma- 
nipulation of its parts. Ideally, the data should flow along 
the bus without any external control [3]. However, with- 
out any dynamic control perfect information transfer for 
most spin chains is possible only if the couplings are pre- 
cisely engineered [H [5J |B] . Although it often suffices to 
engineer the couplings near the end of the chain j7] , fabri- 
cation of chains with such precisely engineered couplings 
is a significant challenge. Perfect state transfer for any 
length of chain is also possible in principle if we have 
full dynamic control over individual couplings between 
adjacent spins, using either dynamic or adiabatic pas- 
sage schemes [9] . However, these schemes require con- 
trol of all couplings, and thus multiple control electrodes, 
which aside from increasing system complexity may also 
be sources of noise contributing to decoherence, an ef- 
fect potentially amplified by the high field strengths and 
longer transfer times required for adiabatic passage. 

Recent work has shown that many systems are con- 
trollable even if we can only modify the dynamics of a 
small part of the system. In particular, the dynamics 
on the first excitation subspace of many spin chains is 
controllable if we can vary the coupling at one end of 
the chain [10], for instance. It has also been shown that 
this global controllability can be exploited to improve the 
transfer fidelities in certain types of spin chains by re- 
peatedly applying certain unitary operations at one end 
of the chain [IT], and in the limit of instantaneous gate 
operations and fast repetition, it can be shown that ar- 
bitrarily high transfer fidelities can be achieved. Theo- 
retical controllability results, however, suggest that one 
could do much better. Almost any local perturbation of 
the Hamiltonian theoretically suffices to effectively con- 
trol the system, and any objective can be achieved in 
finite time by simply switching this perturbation on and 
off at specific times. In this Letter we present a system- 
atic approach to finding the correct switching times to 



achieve fast high fidelity information transfer for various 
spin chains, including both uniform and disordered spin 
chains. We further demonstrate that effective bang-bang 
switching sequences can be found even in the absence of 
a precise model of the system, using adaptive closed-loop 
experiments. This is important as the Hamiltonians for 
a particular physical realization of a spin chain quantum 
wire are usually at best approximately known and sub- 
ject to variation due to fabrication tolerances, although 
this problem can potentially be overcome by experimen- 
tal system identification [12] . 

For proof-of-principle simulations we consider spin 
chains of N spin- 1 particles with a coupling Hamiltonian 

H I = 2 J mn a m a n + J n\n a ni <T l + 4"m' 7 n. ( l ) 

n,rn 

where cr* for * € {x, y, z} are the usual Pauli matrices for 
the nth spin and J mn — J nm are the coupling constants. 
This model covers common types of chains, including XY, 
Heisenbcrg, and dipole-coupled spin chains. For applica- 
tions such as information transfer, it is desirable to re- 
strict the dynamics by ensuring that Hj commutes with 
the total spin operator S = S n =i °n> an< ^ thus conser- 
vation of the total number of excitations. This condition 
is equivalent to xy-isotropy, i.e., J^ nn = J mnl and en- 
sures that the Hamiltonian decomposes into excitation 
subspaces. As usual in work on spin chains for informa- 
tion transfer we assume J mn = J mn and work in the first 
excitation subspace E\. We have dim_Ei = N and the 
off-diagonal elements of the Hamiltonian restricted to E\ 

are Hmn = J mn = J mn f° r m 7^ n - The diagonal ele- 
ments depend on J mn - Most simulations will be based 
on isotropic XY or Heisenberg chains, where J mn = J mn 
and J mn = in the former case, and J mn = J^„ = J 7 mn 
in the latter. For chains with nearest-neighbor coupling 
(nnc) J mn = except for n = m ± 1. An nnc chain is 
uniformly coupled (unnc) if J n .n+i = J f° r an n i an d 
choosing time in units of J" 1 , we may assume J = 1 . 

The objective of optimizing information transmission 
through a spin chain quantum wire is to achieve high- 
fidelity transfer in a minimum amount of time. We aim 
to achieve this by optimizing the switching time sequence 
t = (ti, . . . , tj() of a simple binary switch actuator, such 
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FIG. 1: Population of 10th spin for length-10 disordered nnc 
XYZ chain without control (red) and with optimized bang- 
bang control (blue). Although the actuator only effects a local 
perturbation of the coupling between spins 1 and 2, it signif- 
icantly changes the evolution of the entire system, resulting 
in "constructive interference" of the populations at the Nth 
spin at the target time T = 95.4740. 



as a local control electrode. The evolution of the system 
subject to this bang-bang control is given by 

U(t) = £/«(*!) C/ (2) (* 2 ) ■ ■ • C/ (1) (^-i)C/ (2) (^), (2) 

where U^ m '(tk) = exp(—itkH m ). U^'itk) corresponds 
to free evolution under Hi = Hi (actuator off), and 
U ^(ifc) to evolution under the perturbed Hamiltonian 
i?2 = Hi + He for tk time units. The idea is that al- 
though the actuator changes the Hamiltonian only lo- 
cally in a fixed way, by switching this perturbation on 
and off at different times, we can realize basically in- 
finitely many effective Hamiltonians that give rise to dif- 
ferent evolutions of the system. Most of these will not 
produce a desirable evolution, but some are likely have 
desirable characteristics such as a near unit-height peak 
in the population of a particular state after a short time 
as shown in Fig. [T] The optimization tries to find such 
desirable cases by systematically exploring the parameter 
space, in our case the possible switching time sequences 
t. Specifically, we seek t that minimize the average trans- 
mission error for a quantum bit propagating through the 
chain [3], E(t) = 1 - |(7V|[/(t)|l)| 2 , subject to the con- 
straint that the transmission time T = Ylk=i^k and num- 
ber of switches remain below thresholds T max and K max . 

If the chain Hamiltonian Hi and the perturbation He 
induced by the actuator are known then we can easily 
solve the optimization problem. For fixed K the gradient 
of the objective function is VE = (diE, . . . , 8k E) , where 
the partial derivatives are 



d k E=§^ = -21m [<AW«(t)|l)(l|[/(t)t|AO] (3) 

with f7 (fc >(t) = Ui- ■ ■ H k U k - ■ -U K , where U t is the £th 
factor in ^ . Equipped with this gradient information we 
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FIG. 2: Transfer times and fidelities for Heisenberg spin chain 
with uniform nearest neighbor coupling. Best fidelities achiev- 
able according to in at most 4000 time units without control 
(blue bars / circles) and with optimized bang-bang control 
(red diamonds/line). 



can use either a gradient descent algorithm to find t opt , 
or we can similarly calculate the Hessian matrix of 2nd 
derivatives H = (H k e) with H k i = dkd(,E(t), and use the 
more efficient Newton method with an adjustable step 
size 7 > to itcratively update t 



t s - 7 [i/(£;(t s ))]- 1 vi;(t s ) 



(4) 



until convergence is achieved. In practice there are some 
additional complications as we have constraints on the 
switching times (e.g. t k > i min , £) fc ifc < ^max, etc) but 
these can easily be incorporated in the algorithm [13J. 
The resulting t opt depends on the initial guess t , and 
as with virtually all optimization algorithms, we can 
only guarantee convergence to a local optimum, but the 
method is generally both effective and very efficient. 

We first tested if we are indeed able to achieve fast, 
high-fidelity information transfer using a simple binary 
switch actuator for unnc XY and XYZ chains. For such 
chains perfect information transfer without control or en- 
gineered couplings is possible only for N = 3, and as 
Fig. [2] for the Heisenberg chain shows, the maximum 
transfer fidelity in a limited amount of time decreases 
quickly, and the peak transfer times vary erratically with 
the chain length N. Assuming a simple local actuator 
that switches off the coupling between the first two spins, 
we were able to find switching sequences t that achieved 
transfer fidelities > 0.9999 for both XY and XYZ chains 
up to length 50 with transfer times T(N) « 10N and 
at most K = 4N switches. Even when the precision 
of the switching times t k was limited to 10 -4 , i.e., four 
decimal digits, the error remained below threshold ex- 
cept for XYZ chains with N = 48, 50, where the error 
was slightly above threshold (1.7 x 10~ 4 and 5.6 x 10~ 4 ). 
Surprisingly, we were able to achieve transfer fidelities 
> 0.9999 in many cases even when the switching times 
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FIG. 3: Error log 10 (B(t op t)) as a function of the number of 
switches and target transfer times To, where E(t op t) is the 
minimum error over 10 runs with random initial to 

were limited to only three decimal digits accuracy, and 
for shorter chains even fewer digits appeared sufficient in 
many cases, though at the expense of somewhat increased 
transfer times. This is significant as it shows that we do 
not require "double precision" control of the switching 
times, which would be experimentally infcasible. 

Although the optimal transfer times of T w 107V 
are impressive compared with those for the uncontrolled 
chain, especially considering the transfer fidelities, the 
ease with which the algorithm appeared to be able to find 
solutions suggested that both K and T could be further 
reduced. To investigate this, we systematically varied 
the number of switches and approximate target transfer 
times To for a benchmark case of an XYZ unnc chain of 
length 10. The contour plot of the error on a logarithmic 
scale (log 10 E(t opt )) as a function of K and T (Fig. [3]) 
suggests that for this system near-perfect transfer could 
be accomplished in as little as 50 time units (measured 
in units of J -1 ) with as few as 22 — 24 switches. Push- 
ing the limits tends to slow down the optimization, in- 
creasing the number of iterations and possibly requiring 
several runs with different initial guesses to- Yet, for the 
simple benchmark problem, all runs completed in a few 
seconds (average 1.3s, max. 4.6s) on a standard laptop. 

While chains with uniform nearest neighbor coupling 
provide nice models, in practice the coupling constants 
Jmn are likely to be subject to variation, and direct cou- 
pling may not be limited to nearest neighbors. It is 
thus crucial to investigate the performance of the bang- 
bang control scheme and optimization algorithm for per- 
turbed chains. As Fig. [4] shows, with very few excep- 
tions the algorithm had no difficulty in finding solutions 
that achieved threshold transfer fidelities of 99.99% with 
transfer times of T « 100 and K < 40 switches for mod- 
erately disordered nnc Heisenberg chains of length 10. 
Similar results were achieved for other chains, e.g. XY 
chains. Preliminary simulations suggest that substantial 



FIG. 4: Transfer times T and fidelities T for 100 (nnc) Heisen- 
berg chains of length 10 with randomly perturbed couplings 
Jn,n+i = 1 + e£ (£ Gaussian random variable, (f) = 0, 
CT (£) — 1) f° r 1% an d 10% disorder. Without control, the 
maximum transfer fidelities achievable in < 4000 time units 
and transfer times T are scattered with (JF) = 83.74%, 
(T) = 2, 120 for e = 0.01 (blue diamonds) and {T) = 82.21%, 
(T) = 1, 956 for e = 0.1 (green circles). With optimized bang- 
bang control, near-perfect transfer {T) > 0.9999 for e = 0.01 
and (JF) > 0.9995 for e — 0.1 can be achieved for all 100 test 
systems in approximately 100 time units (red lines). 



gains in fidelities and transfer times are possible even for 
highly disordered chains, although in some cases finding 
solutions becomes challenging and target transfer times 
and/or the number of switches may need to be increased. 

These results are promising in that the simulations 
suggest substantial improvements of both transfer fideli- 
ties and transfer times are possible with very limited 
control and very simple actuators, even for disordered 
chains. However, our optimization procedures assumed 
knowledge of the chain's Hamiltonian and the pertur- 
bation induced by the actuator. A valid objection to 
the practical feasibility of the approach is that this in- 
formation is often simply not available for a particular 
physical system. Furthermore, unknown environmental 
effects may perturb the evolution. Hence, it is crucial to 
consider if we can determine a set of optimal switching 
times for a given system without recourse to an ideal- 
ized model, using adaptive closed-loop experiments. In 
this case, the figure of merit, in our case the transfer 
fidelity, for a particular switching sequence is evaluated 
experimentally as follows: initialize the system, create 
an excitation at one end, apply the control sequence, the 
measure of spin at the other end of the chain. Assuming a 
simple binary-outcome projective measurement, this ex- 
periment is repeated until we have accumulated sufficient 
data to estimate the transfer fidelity to a desired num- 
ber of significant digits. In this setting the optimization 
must find an optimal time sequence based only on the 
limited-precision fidelity measurement data. 
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TABLE I: Performance of various optimization algorithms in 
closed-loop experiment simulations for benchmark problem. 



To assess if we can still find effective switching pulse 
sequences we simulated this situation, i.e., we choose 
a model system to calculate the fidelities, but instead 
of providing the optimization routine with information 
about the actual model as before, it now only had ac- 
cess to an estimate for the fidelity with D significant 
digits for each time sequence t. For our benchmark prob- 
lem (a unnc XYZ chain of length 10) we compared for 
four types of algorithms: genetic, pattern search, Nelder- 
Mead simplex and Newton iteration with discrete gradi- 
ents derived from the limited-precision (simulated) mea- 
surement data. The results are shown in Table |U Pat- 
tern search was abandoned due to extremely poor perfor- 
mance and convergence issues. Most notably, the stan- 
dard genetic algorithm, despite being a popular choice 
for closed- loop optimization experiments in some areas, 
failed completely for this problem. Out of 100 trials 
with different initial populations (50 individuals), not 
one came close to reaching threshold fidelity, and the 
number of fidelity evaluations (#J-), and hence exper- 
iments required, was huge. The simplex algorithm did 
significantly better in that 75% of the trials succeeded in 
finding solutions above threshold fidelity. Application of 
the quasi-Newton-type optimization procedure outlined 
above, requires that the analytical gradients be replaced 
by discrete gradients calculated from the fidelity mea- 
surements. Although the limited precision is a challenge 
here, with careful tuning of the parameters in the discrete 
gradient estimation, the quasi-Newton routine (Newton 1 ) 
far outperformed all other choices in terms of success 
rate, number of fidelity evaluations, execution times, and 
best transfer times. For comparison we have included re- 
sults for our model-based Newton iteration (Newton 2 ) 
using (partially) analytic gradients. As expected, the 
model-based iteration is more efficient in terms of func- 
tion evaluations and execution times, but there is no dif- 
ference in the success rates or average transfer times. 

As estimating the fidelities to a large number of signif- 
icant digits is costly in that it requires many experiment 
repetitions for each fidelity evaluation, we further inves- 
tigated the effect of limiting the precision of the fidelity 
measurements to a few significant (decimal) digits. We 
note that 4 digits is the minimum accuracy required if we 
wish to achieve transfer fidelities > 99.99%. We found 
that with suitable choice of the algorithmic parameters, 
especially with regard to the discrete gradient estimation, 
the quasi-Newton-type algorithm generally continued to 



find solutions with error probabilities < 10~ 4 in almost 
all cases even when the accuracy of the fidelities was lim- 
ited to 5 to 8 digits, and in some cases the algorithm was 
able to find time sequences achieving transfer fidelities 
> 0.9999 for our benchmark problem even when both the 
accuracy of the fidelity measurements and the switching 
times tk were limited four decimal digits. This suggests 
that threshold fidelities can be achieved even with limited 
precision measurements and control. 

We have shown that minimal control using a binary 
switch actuator that induces a local perturbation to a 
fixed Hamiltonian holds considerable promise to improve 
information transmission through spin chain quantum 
wires. The effectiveness of the technique for disordered 
chains and the possibility of model-free optimization us- 
ing closed-loop experiments with limited precision mea- 
surements further enhance its potential appeal. Further- 
more, the technique is not limited to information transfer 
in spin chains with nearest neighbor coupling. Prelimi- 
nary simulations suggest that the technique is still ef- 
fective for more complex spin networks with non-nearest 
neighbor couplings and potentially many other types of 
systems, and that it can be applied to implement control 
objectives other than state transfer, including non-local 
unitary operations. 
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